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ABSTRACT 

Vie develop a theory of Orbits for the Inverse-square central force law which 
differs considerably from the usual deductive approach- In particular, we make 
no explicit use of calculus. By beginning with qualitative aspects of solutions, 
we are led to a number of geometrically realizable physical invariants of the 
orbits. Consequently most of our theorems rely only on simple geometrical 
relationships. Despite its simplicity* our planetary geometry is powerful 
enough to treat a wide range of perturbations with relative ease. Furthermore, 
without introducing any more machinery, we obtain full quantitative results. 

The paper concludes with suggestions for further research into the geometry 
of planetary orbits* 


♦This is a preprint of a paper to appear in he American Journal of Physics . 
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L 5umc QualifyhTC E£ 5 u 


]. l-niraducttCM 

From junior high school on, students of Kknct are caught that Kepler** Laws 
describe the mot inn of planets aroon-d the sun. They are given no hint of how they 
themselves can understand the why of these laws. By high school the students have been 
taught Newton's discovery, that the inverse-square Force law accounts for those beautiful 
ellipses, but rhe connection is not yet for their eyes. After a year or so of college ic’* finally 
urne to plow through the thoroughly stsndardited and unmotivated proofs using intricate 
manipulations with differential equations* 

In this paper we outline an approach to orbital mechanics which i* accessible to 
beginning physics students and presupposes no knowledge of calculus. We give an 
elementary (yet mathematically correct) treatment of Kepler's Laws and also investigate a 
simple first-order perturbation theory for orbits in an inverse-square field. Our theorems 
and proofs arise naturally from trying to understand orbits in terms oF their physical 
invariants, We therefore feet that our treatment provides a better view of "what doing 
physics IS really like than does the Standard route Via algebraic manipulations. 

The key to the method lies in consdeiing the velocity Space picture Tor 4 planet's 
motion about the sun. The concept of a velocity spate is not normally encountered by the 
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itudcnt until he or she is presented with phase space in the context of the formalities of 
Hamiltonlall mechanics. We think there ii little reason for this delay because, though it is 
usually considered an "advanced concept", the velocity space picture of a particle's motion 
lies at the core of Newtonian mechanics, indeed, the qualitative content of Newton's F-ma 
{for our put poses. mAv*FAt) is simply 1 bat physical interactions between objects take place 
by a modification of velocity,, rather than by a change in position, Appreciation of this 
fait Can greatly assist the development Of physical intuition and understanding, and Velocity 
space IS a natural tool for exploring Newton's primarily conceptual breakthrough. In Part 
III, when we discuss perturbations, we thill see what rich dividends can be yielded by 
Soaking at physical phenomena In the right conceptual ftiime-in this case, velocity space- 

i 

In the following presentation, we have tried to walk rather a nartow path between 
two extremes, On one hand, a description of our methods and results would take no more* 
than a few pages if we used the full precision of mathematical apparatus '{including 
calculus) available to science students after a few years of university education. On the 
other hard, we could have spent considerably more spate developing a complete and self- 
contained course for very early physics students. Since wf feet the material can, be useful at 
both levels of physics education, we have attempted a compromise. We apologize both to 
those who find our presentation extended and perhaps verbove, and to those who might 
find It sketchy and incomplete. 


We gratefully acknowledge [he inspiration arid encouragement of Seymour 
PapfrL He introduced us to this way oF thinking about Orb-MS-. and pointed out the basic 
resulti described in Sections t and 3. These sections densely follow parts oF Jiis paper 
tjrefererrce i> which presents a broader yiew of the conception of education in science and 
mathematics from which this work grEw, We would like to thank the editor and referee 
From Tht American Journal of Phjsus for many encouraging and helpful comments and 
also to thank Suiin Jabari of the H.LT Artificial intelligence La hu.ra.tocy for preparing the 
Illustrations for this paper. 

2. The Orbti IS Glitied 

Standard approaches begin with'the arduous task of proving that planetary 
orbits are precise ellipses. We begin by proving a more qualitative proposition that orbits 
are dosed. In doing so w-e dispense with a great deal of analytic clutter, and the important 
special nature of inverse-square orbits which mskes (hem closed Comes into central foots, 

V/c will prove that no orbit like that in Figure I is possible. 

If a planet t rc-jj/j a half-line from the jmt tmia, then a cwiej ft ar 
tAr tome point each Hmf-TWt further Out ar dusrr tn. 


We assume twp pieces of knowledge 


EARTH 



Figure 2; Opposite pieces of orbit, at radii ^ and , 
\ arid in A t] sTid Atg* 


sweep out arEas 




L The force on Lhc planet,, when it is distance r from the sun, is K/'r^ towards ihe 


sun. 

2 - Angular momentum ii Conserved. We use this in the form of Kepler's Law 
that the radius from the sun to a planet sweep* out equal areas in equal limes, This can be 
easily derived and we remind readers of itj simple geometric proof in the appendI k^ 


Mow consider diametrically opposite pieces of the orbit which subtend the same 
(small) angle measured from the sun, as lit Figure 2 . Kepler’s Uw states that 


it 


At- 


What else do we know about the area or r^c time? Geometry Jells us that 


5 

Those r^s. arc too suggestive for US not 10 make a connection with 
F K/r?. In fact 


hence 


Al, 


At 2 3 


and w*r conclude FjAt| ■ F^A;^- Since Fj and Fr> puil in opposite directions, 











We can identify these terms: Awarding to Newton's Second Law, on Each piece of the 
orbit, fAe 11 precisely Av the change in velocity 5 which we cal] the "kick* associated to 
that piere. Thus Che last equation says that the change in velocity over one piece of orbit 
exactly cancels the Change an the opposite piece. Starting on the half-line which the orbit 
crosses twice, divide the orbit ail the way around into similar p^irs of opposite pieces. The 
total change m the planet’s velocity between successive crossings of the half-line is rhe sum 
of the changes in velocity over each small piece; adding these up in opposing pairs, we see 

lha; tile total cnange is tern. Whenever the plan el crosses a given half-line, it has the same 
velocity. 

Now Keplers dictum of equal area in equal time allows us to conclude chat at two 
crossings of the half-line, not only velocity-hut distance from the sun it the same. Figure 3 
shows the areas swept out by the planet m some short time At after successive crossings of 
the half-lme through A, E, and O. The velocities at A and E are equal. Therefor? the 
pieces of cubit AC and ED are boLh equal tot At, but the area of AOC most equat the 
area of EOD, Then A equals B and the orbit closes 

■u 

3. A Thioie si fj: Velocity Spact 

The preceding proof rested on tire fact that the kicks over opposite pieces of the 
orbit have equal magnitudes; 

V’t, = F 2 4t 2 



Figure 3 j ti C"BD t area AOOarea SQD* therefore QA=0B and A-6 



Figure 4: Part of an orbft divided into equal-angle slices. 











equat^n it isn't necessary that the pieces «r Dr bjt be opposite. We Deed 


liut, to derive this 
only that 



and this tl [TUB fa-fry »„ ptew of [ht orb,, ovw which [he tidiif ingle chu ng „ by the 

same small amount. So, if w divide the Whole orbit into small pittas subtending the same 

angle AS , the "kith" veelors for the various pieces all hive the jamt length [figure 4). Not 

only are alt the lengths equal, hot [he totatloo betwe™ suctitive kick vectors is cemstarc and 
equals AS . " 

Thus we have a very simple algorithm for gewriUng Eh* changing velocity as 
the planet moves along its orbit, Starting at a given velocity vector m add on kick vectors 

one after another. Each addition is a step of constant length and successive steps differ in 
direction by the constant turn, A*_ Jr it eas y t9 ^ lhat the algnr|thm . Go FORWARD ^ 

short distance, TURN through a small angle, GO FORWARD the same short distance, 

TURN through the same small angtaU will generate a circle. 5 We concede that the kick 
vectors line up along a circle (Figure 5J. 

We C4h interpret "adding on successive kick vectors" by introducing the notion of 
w/ftrify JpOtt, The velocity of an object is usually described by a “velocity vector* that is, a 




Figure I'iacing the kick vectors end to end. 



direction and a >f«d (length). In comparing different velocities It is useful to put the tail 
of all velocity vectors down at some common point* ft, and to depict a velocity by the point 
where the tip of the velocity vedor lands. With this convention* we on draw two different 
pictures to describe the motion of an object 

_(J) the collection of successive positions of the object In “rear space; and 

{ 2 ) the edition of successive positions of the tip of the object's velocity vector. 
This ls a path in "velocity Space", a picture of how velocity changes. 

i he second picture is called the "velocity space path" or "velocity diagram* 
Figure 6 exemplifies these two kinds of diagrams. 

Velocity Is the thing that changes position* Jitr*j are the things that change 
velocity. To get from an object^ position at one instant, t, to its position at t-At, we add cm 
the vector vAt, To get from the object’s velocity at one instant to the velocity at a slightly 
later time* we add on a kick vector, FAr. Adding up successive vAt vectors gives the 

position spare path: adding kick vectors gives the velocity space path* We can now rejtate 
our result as 

Clrdt Thrtirtm 

For a.0 Cifref rr.avmg in an fitld, tke vtlxiij ffnet 

patch ites en a ctrdt- 



Dn space and velocity spate paths. 


figure 6: One object's pas Hi 



To avoid confusion we point cut that rhc center of this circle is, jwt necessarily at 
the origin in velocity space. 


4- The' Vfloaty Space Path 

Our Circle Theorem tells us that the velocity Space path lies nn a circle. But is it 
a complete circle or just part of it? We can answer that question and a bit more, 

We showed that an orbit which dues manage to get all the way around the sun 
crosses every half-line from the sun exactly once. Such an orbit Is a simple closed curve, 
in a complete revolution, therefore,, the direction of the planet's velocity vector must change 
through 3G0*. (SEP reference -fj That means that fa velactty space aha the path meets every 
half-line from the velocity space origin It follows that closed orbits in position space 
correspond to complete cireks in velocity space, and we have learned, besides, that the 
origin of velocity space Js inside the circle. 

We now have a good qualitative picture of the velocity space path for a dosed 
orbiL (Open orbits are discussed in section IS.} 3n Part IT, we will extract information about 
the position space orbit from our velocity space diagram. 


D 



Figure 7: p = radius x A $ 







H. fnvurtants uI Ihe Orbil 


5 , rfriguitr Mnwtnt&m 

As we have sten, ihc velocity spate path of a planer in an inverse-square field Les 
Dti a circle. One obvious invariant of a eirck ii iis radius. How can we interpret this 
invariant physically? 

W? go* the Cirde in action $ as the result of the algorithm “forward distance 0, 
turn angle Aft, repeat." As one can see (Figure 7), rim generates a circle of radius D/Aft, Ttt 
□ur case Aft was an arbitrary small angle arvd D was the magnitude qf the kick Fit over 
the corresponding small piece of orbit Letting u denote the radius of the velocity Circle,, we 
have 


FA t 

u — —■ 


I* is not immediately obvious that this is a constant. However, wt can Simplify the 
expression using the fact from gEometry that the area swept out over a smalt piece of orbit 
is 


A - — r^A I?. 
2 



A© ~*A 

r 1 by using F*-K/r^ we obtain 







The Term which tells how fast the plan-si is sweeping' out area, is precisely the 

constant called angular momen.rum, L. (See Appendix,} Therefore 

The radius of the velocity circle equals the force constant 
K divided by the angular momentum L: uvKfL 
For a fixed gravitational field, Hie radius of the velocity circle fells us the planet's angular 
momentum: a larger radius gives a smaller angular momentum. 

6 . Onmtatlm 

The velocity circle has another invariant so obvious it is easy to overlook—the 
position of its center in velocity space. As we remarked above, evert though the origin 0 is 
tnsede the velocity circle it need not be at the center of the circle. Let t be the vector 
running from the origin in velocity space to the center of the velocity circle, and let u be 
a radial vector of the circle (Figure S).- Jn terms oft and, ti. we can chink of the planets 
pa.h in velocity space a? follows: at each moment the Velocity "v is the sum of a cojtrfanf 
vector * and a vector u of eomtaru length {equal to u-KfL), V=t*u. The velocity space 
path is generated as the radius ti sweeps around the tip of the invariant vector^ 

There is a <pire remarkable relation between the motion of "r, tbe position space 
radius vector (tail at the sun, head at the planet), and the motion of this "velocity space 
radius* lx! 

C orTtlalim of Angler In Pottmn cm of Vehctty Spate 
At each mmnt the plana j radio? vector r u perpendicular to the 
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Figure The invariant z. 








rtkfriti 11 of tkt vtlMty ClTtle. 

To s-BB this we examine how the kLcks- fit into bath diagrams; In position space each kick is 
parallel to rbe radial vector r. En velocity space the same kick is tangent to the velocity 
circle and hence perpendicular to the velocity circle radius u, 

Hence u is perpendicular tot. 

Correlation of Angle* is a powerful principle, It tells us; 

0) Each point on the velocity spate path corresponds to a unique point in the 
planet's orbit. (The planet cannot attain the same velocity at two different points in its 
orbit.) 

[2) The plant's petition vector sweeps around the sun at the same rate and with 
the same direction (clockwise or counter- clockwise) as its Velocity sw^ps around the Circle 
m velocity space. The two are always 00 degrees out of phase, (Figure 9) 

Now we can.give more meaning to the z vector The planet's speed, v, is 

_S -k 

gfeaEeit when ths u and i vectors are lined up, least when the) 1 are apposed: v mm =u% 
v max“ u+1, Therefore the z VKtoi points in the direction of maximum sperd and 
opposite to the direction of minimum speed. 

The points of greiite&t and least speed occur where the velocity vector is parallel 
to u, and SO perpendicular to the position space radius "r It is not hard 10 show then 
that the point where Speed attains its maximum (respectively, its minimum) corresponds to 
the minimum (respectively,, maximum) distance flORI the Stin. The reader can fill jrr Lh«r 
details, of the proof sketched below: 


-t' 



POSITION 



VELQCJTY 



F^gyre 9: Snapshots of position aftd velocity over art orbit tjf. 














Sfejs 1: At a max imum or minimum distance the velocity v must be 

A 

perpendicular to the radius r. 

Step 2. The velocity diagram shows that there are precisely two points where this 


can occur. 


Step 3. Conservation of Angular Momentum implies that r miJj corresponds to 


V m.n and r mip v mar 


Mow we have some more qualitative information about the shape of the orbit: 
there is precisely one point of maximum distance from the SUrt. and One of minimum 
distance. 1 hey occur on opposite ssdesot the sun since the correspond I ng' 0 Vectors point 
In opposite directions- The i vector determines the orientation el" the orbit. It points in 
the direction of maximum speed, (Figure 10) 


7, $tiapf and SynKidrf 

What more does the length of i tell us about the orbit? Consider what would 
happen if i vanished, t here would then be no direction picked out for maximum speed 
or distance from the SUil. The planet would have to travel around the Sim th a circle at 
unaForm speed. (Another way to see this: v would be etfua! to lT and therefore always 
perpendicular tor and oF constant fcngth'-the characteristic of uniform circular motion.) 

This suggests that z indicates how the orbit deviates from a circle. We can 
make this precise. One obvious measure of the nomciroularity t>F the orbit is the difference 
In the extreme distances from the sun, 


o 



max 


mas 














r maK" r min - 


If we warn an invariant that depends only on the shape and not the size of the orbit it is 
better La se? how much the ratio 


r lftUtf r miri 

diFfers from E. The length aF I, relates the maximum and minimum Speeds: 

W"™ 

To relate speed La distance from the sun we as? angular momentum. If the angular 


momentum is L then 


LwV mn K r ma*“ v rnax * f min 

because at these places in the orbit v is perpendicular to r (Section 6), Therefore 

J^ni gx _ Vm a x „ z 

r min Vmtn ^ u-z 

Since u«K/JU we can also write 

_ K* Lz 
J^tnin Lz 

But K depends only an the nature of the gravitational field so we ace that nur "shape 
invariant" is determined by Li, The larger Li, the more the orbit deviates from a circle. 
In section Jl we will derive the analytic result that the cubit is an dttptf and L i determines 
its eccentricity. 

En fact,, ktiowledge oF K and the velocity diagram essentially determines the 
planct f s motion. The maximum velocity can be read oft the diagram Immediately, as can 
* U-KfLt so we know L. The Shortest radius in position Spate has length Tjnjn^Wv^ and 









u perpendicular to v^ (Section 6). Having chi* one vector,'"r^, we can generate the 
orbit starting at the position tfetetmined by^^ with the following algorithm. 

I Travel a shot! distance vAl in the direction of y, 

2. Measure the change an angle, A#, m position space. 

3. Find the Velocity at this new angle (by rotating "w through Afl and consulting 
the velocity diagram). 

4. Iletorn to step ] This generates the entire-position spacE path. 

Notice that the velocity diagram it symmetric about the line determined by L 
The above algorithm translates this fact into a symmetry of the position space orbit. 
Starting at the nearest point to the tun, construct the orbit in the forward direction for a 
While, along Vj for Afj, then along v 2 for A0 S , and so on. Now go baefc to the 
starting peunt and run the algorithm backwards with the same sequence Df Afl'j, Since the 
velocity diagram ii symmetric we generate the same small segments of orbit, except that they 
have been flipped About the tine perpendicular id t Therefore the entire orbit is 
symmetric abouit this tine. 


S. fumnory 

We have so fit obtained the following information from the velocity diagram. 

1. The radius of the velocity circle determines the orbit's angular momentum: 

u-KfL 

2. The center of the velocity circle determines the orientation of Che Orbit (a 


points in the direction of maximum speed) and its -shape' <Lx determines 

5. From the velocity diagram, we c^n algorithmically rierermihe the whole orbit. 


UK Perturbations 


9. The Firturboilon Fomnia; Rattiai TAraiC 

U ii in the study of perturbiticms, or how orbits change under small kicks other 
than those given by the sun, chat our use pf velocity diagrams really pays off. It pays off 
for a very good reason, which we mentioned in the introduction as a qualitative form of 
Newton's Second Law of Motion; 

Force acts on the paths of particles by changing vtlKUy and not 
poitittm. 

If fail to tike account of ibis fact we may be faced with situations chat appear 
counter'intuLtive. For example, suppose a spaceship in a circular orbit around a planer 
applies a small outward thrust (Figure llajr"' 

How w|H the orbit change? "Intuition" may suggest that the orbit will elongate in the 
direction of rhe thrust, something like Figure life. In fad, the orbit will elongate, but in a 
direction perpendicular to Ihe kick as in Figure 11c, 

To understand this we consider h&w the kick changes the velocity diagram. The 
spaceship started in a circular orbit whose velocity diagram is centered on the origin. 

Since force affects velocity and not position, it is reasonable (and we shall show below) that 
the effect of the kick in velocity space tiaily u just to move the velocity circle in the 
direct boo of the kick (Figure I?). The corresponding change in the position space orbit is 
the 'counter-intuitive - effect described above. 

Our strategy for itudying perturbation; vj|] be to see how kicks change the 



Figure TV: Will the outward kick on orbit {a} produce (b> or (c)7 




POSITION BEFORE 


POSITION AFTER 



VELOCITY BEFORE 


VELOCITY AFTER 


Figure 12: The perturbation 


induced by an outward kick. 

















velocity diagram. More precihHj, we know that th f shape and oriental ion of the orbit is 
determined by LX the t vector times l he angular momentum, SO We want to find the change 
in l£ A(Lz), produced by an arbitrary kick. 

The basic velocity space equation, v -1 . u, fives it - Lv Lu. Since u has 
length KfL we have Lu - Ks p where t is a vector of unit length whose direction is 
determined solely by the object’s position. (Is is perpendicular to the radial vector "t:) To 
compute the effect of a kick At on L* - Lv - K t notice that since Ucki do not affect 
position, 1 is unchanged. K is also unchanged. Therefore the changed Lt is the same as 
the change in Lv, and the first-order approximation to the change in a product of 
changing quantities gives 

Pcrtuth&im Formula: Ad^J-^AL* lAv 

VVe can use the Perturbation Formula to tidy up aur discussion of the "radial 
thrust problem" (Figure li\ Since Av ,j in the radial direction, the angular momentum 
does not change (AL - 0), so the formula implies A(Ll) - LAv. This means that the 
velocity diagram changes from at vector of zero to at vector in the direction of At 
(Figure 12), 

Intuitively, the velocity circle is "pushed" in the direction of the kick. Note that an inward 
kict at Ihe bottom of the position Space Orbit would have the same effect- 



position before position after 




Figure 13i The perturbation induced by 4 tangential kick. 




















l&TaaftitttalTArutii Solar Wfnrf; j^ 2 ’^ l am 


In this section we apply the Perturbation Formula to some other orbit problems. 

Taitgnttat Thnut- Suppose again that a rocket starts to a circular orbit, but this 
tune provides a ling*mi a | kick figure I3>. To determine A(Lt)- vAL * LAv we note that 
vAL is an impulse m the direction of At since Al H positive and v is parallel to Ay, 
Then the neivfy created z vector must he in the direction of the impulse (Figure IS)l The 
c!onpatten in position space is again perpendicular to the kkt 

Tht Solar WtndP Assume the racket IS affected not Only by the planet's gravity 
but also by a small constant force (LtiehrmantVj 'solar wind”), If the perturbing force is 
small compared to gravity, each revolution of the rocket will be nearly ar ellipse. We can 
therefore think of the orbit as an ellipse which varies through timt To compute how the 
ellipse changes we view the wind as providing impulses all along the orbit (Figure H) and 
sum A(Lc) -i vAl + LAtover one revolution. 

The LAv contribution is a net change in the direction of A* To compute vAl 
we notice that AL is positive on the bottom half of the orbit and negative on the top half 
as shown in Figure 35b. Wc can sum the vAL^S by exploiting the symmetry of the orbit 
The vertical components of the vAL's on the [eft cancel the vertical components oF the 
vALs on the right, leaving only a horiwntii component In the direction of (Figure 
IScL This adds with lAv to produce a t vector in the direction of the wind, intuitively, 
the velocity circle gets 'blown" in the direction of the wind, The orbit elongates 
perpendicular to the wind as in Figure 1£, 



Figure 14: The solar wind. 
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Figure 15; The vectors v fa) and t&L (h) for the solar wind. In fc) 
we see that the vertical components of "v^L on the left 
cancel the vertical coh^oti e n ts on the right. 


















Figliire 16: Position space change in the orbit imtf&r the solar-ttfmL 






Since th.E wind does hot afreet the symmetry of the orbit abou? the vertical axis we can 
apply the sime analysis as above to show that the Li vector continues increasing Ifl the 
direction oF the wind. From this we conclude that the orbit becomes more and more 
eccentric while the direction of "a remains constant, and the orbit Continues to elongate 
f>trf>tndictiiar to the wind. 

The orbit becomes closer and closer to a straight lina. and it eventually reaches, a 
point Where the "small" wind can have large qualitative effects over a timescale of less than 
one revolution (the orbit in fact reverses direction), and our method of averting over an 
entire revolution becomes inappropriate. 

The t PaT(t Fiild: IF a is a small constant (we will take it to be positive), we 
can treat the central force field of magnitude r"^ +E ^ as a perturbation of Che r"® field. The 
perturbing force is some force (positive orftegativc) in the radial direction. To understand 
how this perturbation affects the orbit,, we make the important observation That the shape 
of a r ^ or r ^ +t ‘ ! orbit does not depend on (he scale which we use tq measure radius., 
Therefore we can determine shape by using any scale which makes ic convenient to 
compute the effect of the perturbation. For the orbit shown in Figure 17, we scale to make 
the distance OP equal to one. Since 1/r 2 < lfr^ + ^ for r c E and 1/c 2 >■ t/r^ 4 ^' for r > ! the 
perturbing force is as shown. 

For this perturbation the kicks are radial, so L is constant. This means A(Lfc - 
Lit but from the perturbation formula, Al - Q implies AfLl) - LAl?. Hence A t* Av" r 
Now we can sum Aa -Av over an entire orbit The left-right symmetry of the ot bit and 
perturbing force means that the sum of horizontal components of the kicks must cancel. 



Figure 17: Orbit with perturbing fores imdleatsd for r 
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Figure IB: &z perpendicL'lsr to ~z. 
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Figure 19 i Precession of r 


orhi t, 







and thE net At a£ downward, that is, perpendicular to [he original "f Subsequent Aiht Wilt be 
perpendicular to the current "f and this results primarily in a rotation, of A not a charge Jn 
length {Figure 3-3). The "'major axis" or the position spate orbit, by the consequences we 
derived from Correlation of Angie, must follow this counter-clockwise rotation. Though 
the orbit retains its shape Li, it processes {Figure 19). 

Kfarnttf.Tt Should be remarked that in the preceding [WO examples we looked at 
the F vectors as representing Av fur the perturbation formula. Of course wr should have 
used FAt, but, because of the symmetry involved, the At factor can be ignored in those two 
cases. In more complicated situations, though, this does becorne an issue. For cJtimple, we 
invito the reader to us? the techniques of this section lo treat the perturbation Induced by 
an ablate sun, 

r r * r * a 

TrV* ftnalyttc li exults 

II. Tht Orbit ti rt Cmtc Seclio n 

An objection thst IS sure to occur to some of oar readers goes something tike this; 
"All these intuitive methods are fine, but tf you want useful quantitative information you 
have to return to the standard differential equations you've been trying to get along 


Without" 


OF course there are orbit problems our simple methods s^on't handle-. As far the 
standard results, however. are able id derive the orbital equation directiy from our 
velocity diagram using no more than trigonometry: 

The orbit is described in polar Coordinates by the equation 

, L 

u co$Q 

The proof IS a natural correlation of the basic quantities, ft"v, u. [ and L Using 
the definition of angular momentum. At a point in the orbit when theu 1 vector and the"* 
vector differ in direction by an angle £ we construct the angular momentum triangle (see 
append iff )- 


The area of the triangle O AU in Figure SO is by definition Lf2, If h is the height 
□f the triangle tlffen 



Since U and r are perpendicular, the height or the triangle is given by h - 0 - I co& &. 
Therefore 

L ■* rh ■ r ^u - z cos 0}. 

Here p represents the angle in vtltttty spatr between u and ihe fixed vector 7- Correlation 
of Angles implies that 0 also measures the angle in position space from r Ed a flJted Vector 
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Figure £0: The angular ncmentum triangle OAR, 
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figure 21: 5-a.mp'E orbits in position and velocity stac&» 
















perpendicular to *. Therefore r and 0 are potar ccmrdJnai'es In position spare. 

The above equation describes a conic section, When the origin of velocity space 
pie; Wfriun The Cirde. ei>i and the orbit is an ellipse. When the Origin js outride the Circle, ej 
< i and the orbit is hyperbolic. When the circle passes through the origin, u - i and the 
□chit is parabohc ^Figure 2]), 

Writing the equation in the form 


2 



r 



we ger the standard form- for a conic section and see that / K is ihc eccentricity, and 



[2. CfflMtuaUcm of Energy 

Energy conservation does net arise natural iy using this geometric approach 
although we an obtain the result as a simple application. 


Apply the law of tOSinrs to the velocity diagram in Figure 22 to get 





Figure 22 : v £ * z 2 - + u 2 *, 2ui cq$£. 
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“ 2 u z cos# 


V - 2 + u 


Substituting 7 CC3i? *« - y (from Section 11} we obtain 



£ 2 

I + 


2uL 

r 


and rente (since ti - K/L) 



Since 2 and u are constant f« ike of bit, 2 Is a COfllUnt, the total energy, E- It is 

interesting to note that when the planet crosses the semi-minor axis (t perpendicular to"v) "u t 

vV 

X and Y form a right triangle with hence the kinetic energy f\2 is exactly the 

negative of the total energy there. 


13. Ktphr s Third Law 

We can use the relation of angular momentum in area swept out, 2A-Lf h to 
compute the period oF the planet's revolution. ]n one period the planet sweeps out the 
entire area ®r its elliptical urbic. The ares of an ellipse of semi-major a.xts a and 
eccentricity c is given hy A - For the orbit we have 





a "" 2 ( + 

r ) 

max } 

- - f-L + -L.) 



2 \ U~- Z U — 2 J 

The eccentricity is e - zfu, so J | -e 2 * 


C Lu K 

* a ** 

{u% 

'“■Ik 


II 

uj a 

Then the period Js determined by 



LT = 2 J t - e s = 

2 



U ■Jfi 

or 



T - 2 w a % ^ 

2 r 

B % 

L u 

•fir 



In turns o. quantities appearing' in the velocity diagram we get 

T ^ 2*w 

Cu a - Z 2 )% C~2E)% 

K, Open Orbits ' 

father than treat the bypertoht case m detail, we lave the reader to verify the 

following: 

l For an open orbit, the are of the velocity circle which is actually traversed is the 
part shewn below, bounded by the tangents so the circle through the origin of velocity 
space. Velocity space geometry gives Ihe correlation of energy with limiting velocity, v 

















figure Tne velocity diagram for an open orbit. 









-■igune £4: Deflection angle for an op&rt ortsit► 


- -f2 2 “ U 2 -J 2 E {Figure ^3_), 


2. The deflection angle (qngte between the two asymptote of the hyperbola) cin 

be easily found as a function of energy and angular momemunt: tan — -i -V2 E (figure 

, 2 K 

2A) 


15- ^U.ffFJrtehJ for Farlhtr Rtuttrch 

We have by no mean* exhausted the study of the geometry of orbits, in this, 
paper. The geometry of orbits, particularly the perturbation theory, is a rich source of 
problems, even of mini-research projects of the type described! by Luehrmann 6 , Below we 
make some suggestions for problemsand stud} 1 topics. 

An tnrtTUtttVt paradox: The Galilean transformation requires the relation between velociries 
measured in two different flames of refrenra moving With relative velocity Vg to be"v' “~V 

Sl 

4 v 0 ' This it is a simple matter to move into a frame with relative velocity -It and transform 
thei vMtor for an elliptical mbit to ft Why then does one not observe a circular orbit in 
the new frame? £n particular, what fails in the algorithm of Section 8 which dots generate 
a circle in position space given a circle centered about the origin in velocity space? 

An Atd to Astr&galtim? Suppose we had to pilot j. spaceship in a gravitational field (such as 
Simulated in computer “space war" games). Would a velocity diagram be a useful addition 
to our instrument panel? For exampfe, to change from an elliptical Orbit tp a circular orbit 




we need only consult [he velocity diagram and apply a force to cancel the l ’vector- On the 
other hand, lots oF information ii lacking it we use wi ly the Velocity picture. Intercepting 
another spaceship ii a tricky problem Involving timing, (Although merely matching its 
orbit is easy.) What other instruments should supplement or possibly replace a velocity 
diagram'? 

Crmutry of iftt hamimic oitfllatw; A fundamental geometric property of solutions to the 
Jir^ differential equation is that they have a vector temstartt of motion (^or the maximum 
velocity, or the Runge-Lent vector are all possible choices for this constant). The two 
dimensional Harmonic oscillator with equal mode frequencies has a similar structure, 


Solutions have an obvious axis which IDA) also be assigned a magnitude in any number dF 
ways. Can one develop a useful velocity space geometry for fhtit system? Can one treat 
simple perturbation^ as with orbits.? ^ 

Mart Solar Wind.- A further discussion of the solar wind phenomenon could make use at 
the fact that the force field is conservative and therefore the energy of the orbit Is constants 
This implies that the angular momentum decreases as the orbit becomes more eccentric. 
Thus the velocity circle is not only "blown by the wind" but also the radius becomes larger 
and finally infinite when the orbit degenerates to a line. Show that the changing' velocity 
Circle always pastel through two fnced points in the plane. (Figure 

Our method of averaging over entire orbits is only a first-order perturbation 
theory whereas the formula 


d(Lz) 

d s 


d s 


JL 

d-s 
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rien/re Changing velocity circle under solar wind. 













($ is any para meter) i& exact. One Would flke la hive a treatment at the SO Jar ■wind Which 
Works near fhe turnaround paints, i, r u approaching I. 

FinaTly, is there a complete perturbation theory based pn the geometry of orbits? 
In particular how can one treat perturbations out of the plane of the orbit? 


?!pf)£nd[K; Angular Me*icTilii.r7: ufirf Kepler's Sramd La if 


Throughout this paper we have been assuming Kepler's Second Law. There is » 
simple geometric proof of this which we place here in an appendix because it did not 
originate with its. It can be found in Newton's Prirtcipia! 

The angular momentum, which we denote by L, is defied to be twice the area of 
rhe triangle determined by the velocity vector and the radius vector from the sun to the 
phneL As shown below L is constant if the velocity doesn't change - the triangles have 
equal areas Since they have equal bases (the length of v) and equal heights (Figure At). 
More remarkably, L Ttmaills wnitont If me change the velocity by applying any kick tn the 
radial direction f/emwtfj or may from The affect of a kick Ay on the angular 

momentum triangle is itinerated below (Figure A2). The kick changes v tov' but triangles 
operand OPQj have the same base, OP, and the same height (h in the diagram) since OP 

and Qg are Parallel. Therefore OP&ami OPQj have the saw area, and angular 
momentum Js unchanged, 

A planet moving about the sun, subject ic m force but the sun's gravitation, his 
every kick applied in the radial direction. None of these change the angular momentum. 
Which Is therefore an to variant of the planets orbit. To find a geometric interpretation of 
this fact, we examine the orbit at time intervals it small enough that the Velocity docs not 
change much over each interval. In each interval the radius vector sweeps out a small 
triangle. The area of am of these small triangles is 



Figure A!: Triangles o£ equal area 












Figure A2i Area OFQ » JVrea QFft' 





Ihaal * ^gwht) ■ j. u 1 r a i 

' 7 — v A t k h=~ L A t 
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( J i as in Figure AlJ. The total area. nvtpt out over sonic tong tinte 


T T-itj ■* 

is i he i’urn of the areas of the small triangles 




JL ^ 1 

A = "Lit ^—- LiLt... 

2 1 2 * 

=■- i L (4t i+1 , 2 *...) 

" L T 

2 

This gives Kepler's Second Law 

For a tody moving in a rad Lai force field, the radial vector sweeps 
out equal areas lh equal tim^ 

It is unfortunate that this proof Ls rot more often presented in physics courses (although 
Feynman ^ discusses stand there is a movie demonstrating this argument). 
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